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We describe a method of white-beam inelastic neutron scattering for improved measurement 
efficiency. The method consists of matrix inversion and selective extraction. The former is to resolve 
each incident energy component from the white-beam data, and the latter eliminates contamination 
by elastic components, which produce strong backgrounds that otherwise obfuscate the inelastic 
scattering components. In this method, the optimal experimental condition to obtain high efficiency 
will strongly depend on the specific aim of the individual experiments. 

PACS numbers: 



I. INTRODUCTION 

Inelastic neutron scattering has come to be recognized 
as indispensable in modern materials science, because a 
material's spin and lattice dynamics provides unique in- 
formation about a system's Hamiltonian. A complete 
description of these excitations in momentum (Q) and 
energy (E) space is needed to fully reconstruct the inter- 
actions that govern a material's behavior on the atomic 
scale. However, the technique normally requires a large 
volume of sample, often on the order of several cubic cen- 
timeters 1 -. This is a very significant limitation in research 
to develop new materials with novel functions. 

In neutron diffraction experiments, however, the de- 
velopment of time-of-flight (TOF) technique allowed the 
use of a white beam for increase in measurement effi- 
ciency compared with the conventional method using a 
monochromatic beam. Each wavelength (energy) com- 
ponent can be resolved by TOF, as shown in Fig. HJa), 
and finally merged into a single diffraction |Q| pattern 
for a powder sample (time focusing) or a Q map for a 
single-crystal sample. Unfortunately, the same cannot 
be applied to inelastic scattering, because the different 
Ei (incident energy) and Ef (final energy) components 
are entangled at the same TOF for the same pixel on the 
detector, as shown in Fig. [ljb). Thus, it has been con- 
sidered that either Ei or Ef must be monochromatized 
or must be analyzed, either of which incurs a large loss 
in neutron intensity. 

Another remarkable method, called cross-correlation, 
was developed over four decades ago as an extension to 
the white-beam diffraction 2 '. The method basically in- 
volves extracting the elastic components and removing 
the inelastic components 3 . As shown in Fig. f2J a special 




FIG. 1: (Color online) Plot of TOF against position with- 
out sequence chopper. Solid arrows indicate the most proba- 
ble neutrons. Dotted lines sectionalize each Ei channel, (a) 
Diffraction, (b) Inelastic scattering. 



mechanical chopper modulates a white incident pulsed 
beam with a pseudorandom open/close sequence, and N- 
times cyclic phase shifts of the modulation generate a set 
of N data with different Ei contrast. Then, on the basis 
of the contrast, the data for each Ei can be mathemati- 
cally resolved. 

The mathematical formalization is given below. Here, 
for convenience, parameters and functions are renamed 
and redefined from those in the original papers^ 3 . The 
intensity detected at a specific TOF at a specific pixel of 
the detector, J bs(p) (p = 1> • • • } N), is described by 



-fobs(p) 



3=1 



(1) 



where p is the phase shift of the sequence (phase of se- 
quence chopper), F(k) is the fc-th element in the sequence 
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setup. The proposed method has two novel aspects: the 
introduction of an inverse matrix representation and a 
proposed method for selective extraction. The former af- 
fords a different solution to Eq. (|TJ) to resolve each E\ 
component in the white-beam data. The latter elimi- 
nates contamination by elastic components, which other- 
wise produce strong backgrounds. Finally, we present es- 
timates of some instrumental specifications for the TOF 
polarized neutron spectrometer, POLANO, being con- 
structed at J-PARC. 



FIG. 2: (Color online) Example of TOF-position diagrams 
with sequence chopper (N — 5) at sequence phases p = 1 (a) 
and p = 2 (b). By cyclically shifting the phase, a raw data set 
of /obs(l), Jobs (2), Jobs (3), Jobs (4), and 7 bs(5) is obtained at a 
specific TOF and at a specific pixel on the detector. After the 
measurements, each Ei component of 1(1), 7(2), 7(3), 1(4), 
and 7(5) can be resolved mathematically. 

F consisting of only (close) and I (open), F(k + N) is 
defined to be equal to F(k) for k = I, . . . , TV, j is an index 
for E[, is the intensity coming from the j-th Ei com- 
ponent in a white incident pulsed beam for j = 1, . . . , TV, 
and B is the background. The pseudorandom sequence 
F is restricted by 

TV = 2" - 1 (n: integer), (2a) 
F'(k) = 2F(k) ~ I, (2b) 

JV 

XV(*) = 1. (2c) 

k=l 

N 

F'(k)F'{k + k') = (TV + l)5 , k ' - 1- (2d) 

fc=i 

This type of sequence F is currently called a maximum 
length sequence, which is generated by a simple recur- 
rence formula and is widely applied in the field of digital 
communications 4 . Combining the above equations, one 
can resolve each E\ component, 

= ^En* +p)Job.(p) - j^jb. ( 3) 

It should be noted, however, that the method can- 
not be directly applied to inelastic scattering because 
the elastic components and their large statistical errors 
obfuscate the very weak inelastic components. This is 
probably why the method has not been realized thus far 
in an actual instrument dedicated to inelastic scatter- 
ing. In fact, for the new CORELLI instrument under 
construction at the Spallation Neutron Source (SNS) in 
Oak Ridge, Tennessee, the method will be used mainly 
to study diffusive elastic scattering such as in frustrated 
systems and ionic conductors^. 

This paper presents a modification to this method for 
a more practical white-beam inelastic neutron scattering 



II. INVERSE MATRIX REPRESENTATION 

We formalize an alternative solution to Eq. ([T]) without 
the conditions Eqs. (f2"a|) - (|2dl) . Here, the measurement 
of /obs(p) (p = 1,...,N) is the same as in the original 
method except for the kind of sequence. Ignoring B for 
simplicity, Eq. ([1} can be represented by 

-fobs = FI, (4) 

where JT bs is the vector [7 b s (p)] (p = 1, . . . , TV); F is the 
matrix [F p=1 ,F p= 2, . . . , jF p =jv]; F p is the sequence vec- 
tor [F(j +p)} (j = 1,. . . , N); and I is the vector 
Hence, one can resolve I by 

I = F- 1 W (5) 

Consider, for example, the sequence (0, 1, 1, 0, I) at 
p=l: 
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Hence, one can obtain 



r/(i)- 




" -I 


-1 


-1 


2 


2 " 




" -fobs(l) 


1(2) 


1 

~ 3 


2 


-1 


-1 


-1 


2 




7obs(2) 


7(3) 


2 


2 


-1 


-1 


— I 




7obs(3) 


7(4) 
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Thus, almost all types of sequences can be used as long 
as F^ 1 exists. Taking into account B again, one can also 
identify a sequence to minimize \F~ 1 (B, B, B, B, B)\, for 
example, by trial and error with many numerical trials. 

It should be noted that this general matrix formal- 
ization is not considered superior to the maximum 
length sequence. However, the general matrix formal- 
ization does afford an advantage when the conditions of 
Eqs. (|2"a|) — (|2d[) arc not satisfied on actual instrumenta- 
tion, for example, because of insufficient switching speed 
between and 1 for high E{ range or high resolution. 
In this paper, we use the general matrix formalization 
only because selective extraction, proposed in the next 
section, also does not fulfill the conditions. 
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III. SELECTIVE EXTRACTION 



IV. STATISTICAL EFFICIENCY 



We explain the proposed selective extraction method 
using the above example in Eq. ([5]). First, one needs 
to prepare another chopper with the inverted sequence- 
from open/close to close/open, that is, from 1/0 to 0/1. 
The inverted chopper gives another data set, 



(8) 



where [J bs{p)] is the raw data obtained for the phase 
p. Then, we consider the case where 7(3) is the elastic 
component for the targeted TOF and pixel of the detec- 
tor. Our purpose is to remove 7(3). Thus, by selectively 
extracting only the arrays in which the third column is 
(sequence chopper closed) from Eqs. (JB]) and ©, one 
can reconstruct a good quality data set: 



" Job S (l) " 




■ 1 
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1-/(1)1 


Jobs(2) 
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7(2) 


Jobs(3) 
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Jobs (4) 
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" 1 








1 


0" 




1-/(1)1 


Jobs(2) 







1 








1 




7(2) 


7obs(3) 
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(9) 



Hence, 



Jobs(l) 
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Jobs (2) 
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7 bs (3) 
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Jobs (4) 
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. 7 obs (5) _ 
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/(I) 

7(2) 
7(4) 
7(5) 



(10) 



This equation can be solved by dropping one array and 
using the inverse matrix, or by the least squares method. 
Further, we can reincorporate the background term B: 



(11) 
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Hence, 



7(1)" 
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1 - 




Jobs(l) 
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1 
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7(4) 
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-1 . 
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(12) 



We emphasize that the selective extraction method is 
applicable not only when the third channel is elastic but 
also when an arbitrary channel is elastic. One can remove 
the elastic components at all TOFs and pixels with only 
the two data sets. 



Price and Skold reported that the cross-correlation 
method is not always better than the conventional 
monochromatic method in terms of statistical efficiency 6 . 
The reason is that white-beam data counts [7 obs (p)] 
and [Jobs(p)] are inevitably large and are accompa- 
nied by large statistical errors [A7 obs (p)] = \J 7 obs (p) 



and [A J obs (p)} 



which are propagated to 



y/Jobs(p) 

with further enhancement by adding and subtract- 
ing [7 obs (jj)] and [J ohs (p)}. For example, {A7(l)} 2 = 
(-1/3) 2 • 7 obs (l) + (-1/3) 2 • 7 obs (2) + (-1/3) 2 • 7 obs (3) + 
(2/3) 2 -7 obs (4) + (2/3) 2 -7 obs (5) in Eq. ©, and {A7(l)} 2 = 
(-1) 2 • Jobs(4) + l 2 • 7 obs (5) in Eq. O, where A7(l) 
denotes a statistical error of 7(1). In addition, since N 
measurements with cyclic rotation of sequence are needed 
to resolve [7 obs (p)] to [7(j)] with only one-measurement 
statistics, the use of white beam does not overall mul- 
tiply the measurement efficiency so much. In the use of 
maximum length sequence, a statistical advantage can be 
obtained only for special j(s) (E[ channel(s)), in which 
the signal of interest is more than twice the average 
counts per channel 6 . Hence, the cross-correlation meth- 
ods would be suitable only in the cases of phonon reso- 
nance, magnon resonance, or elastic scattering. 

The above situation is essentially the same as that 
in our modified cross-correlation method with selective 
extraction. Here, we roughly estimate the statistical 
efficiency based on the assumption that consists 
of only inelastic scattering counts with similar magni- 
tude and without huge elastic components. First, in the 
cross-correlation methods without selective extraction in- 
cluding the original one and the aforementioned general- 
matrix one, the error of I^ nsl \j) is described by 



{A/ (ns 1)(?)} 2 = ^ (i ^ )2/ (ns 1)(p) 
r (nsl)v 



(13) 



a j (/obs ') ~ -^opon«j 



</ ( 



nsl)\ 



where the superscript (nsl) denotes the non-use of se- 
lective extraction, A' 1 corresponds to the inverse ma- 

r(nsl) 



trix, a) = Y,p=i(Ajp) 2 , = E^IZ'(P)/N, 



</ ( 



open is the number of opening channels (~ N/2), and 

nsl)\ _ 



J2f =1 I {nsl) (j)/N. Hence, the statistical effi- 



ciency 77^ nsl ' ) is estimated by 



(nsl) 



s/2c 



Nai 



J (nsI) (j) 
A7(" sl )(j) 

■vV M) > 

• v< /(nsl) >< 



(14) 



4 



where I^ nsl \j) = Cj(/( ns1 '). Next, in the modified cross- 
correlation method, the error of I^ si '(j) is described by 



JY 



P =i 



N 



E(^p) 2 {^ 1) (p)/2} 



(15) 



P =i 



^ Ia 2 , /r( nsl )\ _ ijv b 2 (I^) 

2 J '\ obs / 2 °P cn J \ /' 

where the superscript (si) denotes the use of selective 
extraction, K^ s (p) = J<£(p) or jW(p) like in Eq. 

^obib) ~ ^ob B (p) ~ ^obs'^P)/ 2 since selective extraction 
involves the use of two data sets: one from the original 
chopper and the other from the inverted chopper, B^ 1 
corresponds to the final inverse matrix like in Eq. (fT2|) . 
b 2 , = J2 P =i(-Bj> p ) 2 i an d j' denotes the row number in 

B- 1 with which I {sl) {j) can be obtained (e.g., j' = 1 for 
j = 1 and f = 3 for j = 4 in Eq. (fT2"]l). Hence, the 

statistical efficiency r]j is estimated by 

(sl) = = J(™'>(j)/2 



(16) 



2\/ A'open^j' 




• \/< j(nsl) )» 



where I (sl) 0') = / (nsl) (j)/2 since selective extraction in- 
volves the use of two data sets again. Then, in the 
conventional monochromatic experiments, the efficiency 
i)] moll °' required to obtain the same data set of all 
(j = l,...,N) over the same total measurement time 
using monochromatic beam is estimated by 



(mono) 



.(/(nsl) 



^ = Vc--^7M)), (17) 



which corresponds to the measurements described by 
A = E (identity matrix) in the general matrix formal- 
ization. Thus, the ratios of the statistical efficiencies are 
estimated by 



n 



(nsl) 



(mono) 



v 



1 2cj 



V 



(si) 



(mono) 



v 



2 • 



(18a) 
(18b) 



Using example sequences summarized in the Appendix, 
we numerically calculated the values of a,j and bj' and 
the criteria of Cj to retrieve the of interest more 
efficiently than the use of conventional monochromatic 



beam. For the original cross-correlation method with 
maximum length sequence without selective extraction, 
the results are a 2 = 0.44, 0.23, 0.12, 0.062, 0.031, and 
0.016 for all j at N = 7, 15, 31, 63, 127, and 255, re- 
spectively. Hence, (r^ /f^ 0101 "^) > 1 gives the crite- 
ria of c 3 > 1.5, 1.8, 1.9, 2.0, 2.0, and 2.0, which are 
consistent with Price's criterion of Cj > 2 for all the 
TV; the of interest should be more than twice the 
average counts per channel^. For the modified cross- 
correlation method, b 2 , = 0.88, 0.39, 0.22, 0.15, 0.096, 
and 0.059 (f ^ N: I{j)) and b% = 0.28, 0.31, 0.63, 0.78, 
0.20, and 0.30 (j' — N: constant background). Hence, 

(if'/lj™ 01 ) > 1 U ¥= N '- IV)) g^es c* > 6.2, 5.8, 6.7, 
9.4, 12, and 15. The latter criteria are harder than the 
former ones in the absence of huge elastic components, as 
expected. In addition, since Cj increases as iV increases 
in the latter, it would be better to set N less than about 
60. 



In this way, both the original and modified cross- 
correlation methods can efficiently give only the 
components with relatively large Cj among [/(j)]. There- 
fore, it is important to tune the experimental conditions, 
such as the ranges of E- u TOF, and pixel used, such that 
the components of interest become the strongest in in- 
tensity among [J(i)]. In this sense, we would like to 
emphasize that the modified method with selective ex- 
traction can remove not only an elastic channel but also 
an arbitrary one channel of no interest with relatively 
strong intensity among [/(?)]; for example, spurious neu- 
trons scattered on unexpected paths and spurious neu- 
trons coming from the previous flame. 

In practice, (1) for elastic scattering, the original 
method without selective extraction could be safely ap- 
plied, as has been reported in the pastSi2iSi£, since elastic 
scattering is normally the strongest among [/(j)]. How- 
ever, the modified method might improve the data as 
an insurance to remove a spurious channel. (2) For 
quasielastic scattering and low-energy inelastic scatter- 
ing, which will have the next strongest intensity and will 
necessarily overlap with the strongest elastic scattering 
in TOF, the modified method will be effective to remove 
the elastic scattering. (3) For higher energy modes with 
relatively weak intensity, one must remove not only an 
elastic channel but also quasielastic and low-energy in- 
elastic channels, which contaminate the higher energy 
data as a source of statistical errors. Therefore, it would 
be better to narrow the Ei range from white to quasi- 
monochromatic, which will allow us to avoid all the elas- 
tic, quasielastic, and low energy inelastic scatterings by 
TOF, in conjunction with the modified method to remove 
an arbitrary channel of no interest again. In any case, 
to efficiently utilize the cross-correlation methods, users 
should recognize that the optimal experimental condition 
will highly depend on the specific aim of the individual 
experiments. 
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V. SPECIFICATIONS FOR PRACTICAL 
IMPLEMENTATION 

We are constructing a TOF polarized neutron spec- 
trometer called POLANO at a decoupled moderator at 
J-PARC. Because polarization devices such as Heusler 
crystals and spin niters lose a large proportion of neu- 
trons, we need a method to gain an increase of over a 
factor of 10 in measurement efficiency. Our method of 
white-beam inelastic neutron scattering is one candidate, 
though it is applicable to both unpolarized and polarized 
inelastic neutron scattering. In this section, taking this 
spectrometer as an example, we present some specifica- 
tions for the sequence chopper. 

For this system, we set the distance between moderator 
and sample L\ as 17.0 m, the distance between sample 
and detector L 2 as 2.0 m, and the distance between se- 
quence chopper and sample L3 as 2.0 m. The time width 
values at the decoupled moderator at J-PARC, Ai m , were 
used for each E\ (Table HJ. For each Ei, the energy resolu- 
tion AE/Ei was kept < 0.04. Then, using the analytical 
formula for energy resolution^, we evaluated the opening 
time per Ei channel required at the sequence chopper, 
Ai sc i. 

Also, on the assumption that a sequence chopper is 
alternately closed without generating a sequence as an 
example, as shown in Fig. [3J we estimated the condition 
of At sc2 necessary to avoid contamination by the elastic 
tails from neighboring E\ channels. The condition can be 
estimated by Aid < Ai c h, where Aid is the time width of 
the elastic component at the pixel on the detector with 
E\ , and Ai c h is the time difference between two TOFs of 
neighboring E\ channels at the pixel, as defined in Fig. [3] 
The two parameters are described by other parameters, 



At d Atd 



A< d 



Li 



At„ 



At m +At s 



;(Lx-L 3 ) 



At 



At. 



Li 



At„ 

L 2 



\-At 



ch 



(Li 
{2At sc2 ) 



L 3 ) 



At m ,(19a) 



(19b) 



As summarized in Table U there exists a solution for 
Ai sc . For example, over a wide E\ range of 10-80 meV, 
At sc — 9 //sec simultaneously satisfies the constraint of 
AE/Ei < 0.04 and the avoidance of the tails from neigh- 
boring Ei channels. This opening time can be realized 
by a set of counter-rotating disk choppers with the fol- 
lowing parameters: 700 mm-0, 20 mm/channel, and 350 
Ha£. In addition, for the sequence chopper, note that 
an arbitrary sequence can be generated by printing a se- 
quence clockwise and then counterclockwise on counter- 
rotating disk choppers. The sequence chopper generates 
0.5 • (tt • 700 (mm))/20 (mm) = 55 Ei channels. 



VI. SUMMARY 

We developed a modified cross-correlation method for 
an increase in measurement efficiency of inelastic neutron 




sequence 
/I J/1 H I I H chopper 
I (counter- 
rotating) 

(At m +At sc ") (L '" L3) . 
v ' moderator 



tof 



FIG. 3: (Color online) Diagram showing TOF and posi- 
tion with sequence chopper in the estimated specifications for 
POLANO. Solid arrows indicate the most probable neutrons, 
and dotted arrows indicate most inaccurate neutrons. The 
sequence chopper is alternately closed using a set of counter- 
rotating disk choppers. The alternating closing parts do not 
generate the sequence. 



TABLE I: Results of numerical estimation of system specifica- 
tions for POLANO. The values At sc i and At SC 2 are obtained 
so as to satisfy AE/Ei < 0.04 and avoid contamination by 
the elastic tails from neighboring Ei channels, respectively. 
All the time widths are defined as full widths at half maxi- 



Ei (meV) /At m 


(/tsec) At. 


ci (/xsec) 


At sc2 (psec) 


5.0 


50 


< 31 


> 11 


10 


32 


< 22 


> 7 


20 


20 


< 16 


> 4 


40 


13 


< 11 


> 3 


80 


9 


< 9 


> 2 



scattering. First, a different solution using an inverse 
matrix representation was formalized to resolve each Ei 
component in the white-beam data. Second, a method 
of selective extraction was proposed to avoid contamina- 
tion by elastic components. Third, taking spectrometer 
POLANO at J-PARC as an example, practical specifica- 
tions for the sequence chopper were estimated. Experi- 
mental situations giving high efficiency, however, would 
be quite limited and complex. 

The example sequences used for the calculations in 
Sec. 4 are summarized. For the original cross-correlation 
method, (1, 0, 0, 1, 1, 1, 0) at N = 7, (1, 0, 1, 1, 0, 0, 1, 
0, 0, 0, 1, 1, 1, 1, 0) at N = 15, (0, 1, 1, 0, 1, 0, 0, 1, 0, 
0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0) 
-At N = 31, (1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 
0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 

0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 

1, 0) at N = 63, (1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 
1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 
1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 
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o, 


0, 0, 0, 1, 


0, 1, 1, 1, 1, 1, 0, 


0, 


1, 


0. 


1, 


0. 


1. 


1, 1, 0, 


0, 


1. 


1, 0, 1, 0, 


0, 0, 1, 0, 0, 1, 1, 


1, 


1, 


0. 


0. 


0. 


1. 


0, 1, 0, 


0, 


o, 


0, 1, 1, 0, 


0, 0, 0, 0, 1, 0, 0, 


0, 


0, 


0. 


0. 


1, 


1. 


1, 1, 1, 


1, 


1. 


0) at N = 


= 127, (1, 1, 0, 1, 


1. 


0. 


0, 


1. 


1, 


1, 


1, 0, 0, 


0, 


1. 


1, 0, 1, 0, 


1, 1, 1, 0, 0, 1, 0, 


0, 


0, 


0. 


1, 


1. 


1. 


1, 0, 1, 


1. 


1. 


0, 1, 1, 1, 


1, 0, 1, 0, 0, 0, 0, 


0, 


1, 


0. 


0. 


0. 


0. 


0, 0, 1, 


0, 


1. 


1, 0, 1, 1, 


1, 1, 1, 0, 0, 1, 1, 


0, 


1, 


1, 


1, 


0. 


0. 


0, 1, 0, 


1. 


1. 


1, 0, 1, 0, 


0, 1, 1, 0, 0, 1, 0, 


1, 


0, 


1. 


0. 


1, 


0, 


0, 1, 0, 


0, 


1. 


0, 0, 0, 1, 


0, 1, 0, 0, 0, 0, 1, 


0, 


0, 


1. 


1, 


0. 


1, 


0, 0, 0, 


1, 


1. 


1, 1, 1, 0, 


1, 0, 1, 1, 0, 1, 0, 


0, 


1, 


0. 


1, 


0. 


0. 


1, 1, 1, 


1. 


1 

1. 


i n n n 

1, 0, 0, 0, 


0, 0, 0, 0, 1, 1, 0, 


1, 


1, 


0. 


1, 


0. 


1. 


0, 0, 0, 


1. 


o. 


0, 1, 0, 1, 


1, 1, 1, 0, 0, 1, 0, 


1, 


1, 


0. 


0. 


0. 


1. 


0, 0, 0, 


1. 


1. 


0, 0, 1, 1, 


0, 0, 0, 1, 1, 1, 0, 


0, 


0, 


0. 


1, 


1, 


0. 


0, 0, 0, 


0, 


1. 


1, 1, 0, 1, 


1, 0, 0, 0, 0, 1, 0, 


1, 


0, 


1, 


1, 


0. 


0. 


1, 0, 0, 


1. 


1. 


1, 0, 0, 1, 


1, 1, 0, 1, 0, 1, 0, 1, 


1, 1, 1, 


1, 1, 1, 1, 0) at 



N = 255. 

For the modified cross-correlation method, the original 
chopper sequences are (1, 1, 1, 1, 0, 1, 1) at N = 7, (1, 

0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1) at N = 15, (1, 0, 1, 

0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 

0, 0, 0, 1, 0, 1) at N = 31, (1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 

1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 
1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 
0, 0, 0, 0, 1, 1, 1, 1) at N = 63, (0, 1, 1, 0, 1, 0, 0, 1, 0, 
0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 
0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 
0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 
0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 



0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 

0, 0, 0, 0, 1, 0, 0, 0) at at N = 127, (0, 0, 1, 1, 0, 0, 1, 0, 

1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 
1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 
1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 
1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 
0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 
0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 

0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 

1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 
1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 
0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 

0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 

1, 1, 1, 0, 1) at N= 255. 
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